More on Gribov copies and propagators in Landau-gauge Yang-Mills theory 
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Fixing a gauge in the non-perturbative domain of Yang-Mills theory is a non-trivial problem due 
to the presence of Gribov copies. In particular, there are different gauges in the non-perturbative 
regime which all correspond to the same definition of a gauge in the perturbative domain. Gauge- 
dependent correlation functions may differ in these gauges. Two such gauges are the minimal 
and absolute Landau gauge, both corresponding to the perturbative Landau gauge. These, and 
their numerical implementation, are described and presented in detail. Other choices will also be 
discussed. 

This investigation is performed, using numerical lattice gauge theory calculations, by comparing 
the propagators of gluons and ghosts for the minimal Landau gauge and the absolute Landau gauge 
in SU(2) Yang- Mills theory. It is found that the propagators are different in the far infrared and 
even at energy scales of the order of half a GeV. In particular, also the finite-volume effects are 
modified. This is observed in two and three dimensions. Some remarks on the four-dimensional case 
are provided as well. 

PACS numbers: 11.15.Ha 12.38.Aw 14.70.Dj 
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I. INTRODUCTION 



Green's functions encode the full dynamics of a the- 
ory In gauge theories these are in general gauge- 
dependent, and it is necessary to determine them in a 
fixed gauge. In the perturbative regime, this can be con- 
veniently done by local gauge conditions @. However, to 
determine the Green's functions over the full momentum 
range it is necessary to extend the definition of a gauge 
into the non-perturbative domain. Unfortunately, there 
is no unique local possibility to extend the local gauge- 
fixing conditions of perturbation theory into this regime. 
One is therefore faced with a, possibly innumerable, num- 
ber of gauges all corresponding to the same perturbative 
gauge. This will be further discussed in section |TT] and 
IV Al for the example of Landau gauges. In particular, 
two different Landau gauges will be investigated here in 
detail, the minimal Landau gauge and the absolute Lan- 
dau gauge. 

The importance of differing non-perturbative comple- 
tions of perturbative local gauge conditions is that they 
potentially lead to differences in the gauge-dependent 
Green's functions. Of course, these differences have to 
become subleading at sufficiently large momenta, where 
non-perturbative effects become sub-leading. However, 
and as will be shown in section ITVl these differences may 
be qualitative at low momenta, and may still be quanti- 
tatively relevant at intermediate momenta of the order of 
the scale Aym and above. Hence, not only calculations 
in the asymptotic infrared momentum range have to take 



this into account. Precision calculations of observable 
quantities, which are based on gauge-dependent correla- 
tion functions, may also need to take care of this. Apart 
from the practical purpose of calculating the Green's 
functions, this problem is also of a more fundamental 
interest. Confinement scenarios, like the one of Gribov 
and Zwanziger 0, and the one of Kugo and Ojima 
7] (GZKO), give predictions for correlation functions in 
certain, well-defined non-perturbative gauges. To check 
these predictions and hence the scenarios, it is neces- 
sary to determine the Green's functions in exactly these 
gauges. 

The determination of the Green's functions is not a 
trivial problem. To obtain them in the full momentum 
range, various methods have been emplc 
methods, like Dyson-Schwinger equations 
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functional renormalization group equations 
32| , stochastic quantization approaches 0, [f| 
or effective action methods [H, [U, [H, Us H3] have to 
rely on assumptions, in general. Latti ce g auge theory 

MMJMJM lliQiiMliiliEIIElial^l^lHl, 

|52| . l53l . l54l l5a ] is not in need of such assumptions, but 
suffers from artifacts, in particular finite volume effects 
SI B HI E3, HI HI- On top of these problems comes 
the necessity to fix the same non-perturbative definition 
of a gauge in all methods. For functional calculations 
this can possibly be implemented by means of boundary 
conditions 29]. In case of lattice gauge theory, this is 
a hard numeric al p roblem, which has been investigated 
previously HE EE HE HE [EH following the pioneering 
work [5^], and will be discussed in great detail here as 
well. 

Therefore, neither approach alone has been able to pro- 
vide a final result for the Green's functions. For a reli- 
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able comparison therefore all of these problems have to 
be addressed. However, the uncertainties due to the in- 
herent problems in this comparison have sparked quite 
a prolonged debate on the infrared behavior of Green's 
functions, and whether they are in agreement with the 
predictions of the GZKO scenarios or not, for a summary 
see [29[. In particular, the investigation with lattice cal- 
culations has, so far, not yet provided a clear picture. 
This will be discussed in more detail in section [TVl and 
the issue of gauge fixing on the lattice is covered in sec- 
tion [HT] and appendix IB1 



In addition to the problem of gau ge-fi xing, finite vol- 
ume artifacts are the main problem |43l . 0, [4^, H3, [H, 
l5r| . The influence of these artifacts has been investigated 
in detail [H, 0, H^, |4?], HH , and seems to indicate a con- 
tradiction between the GZKO scenario and the actual 
form of the Green's functions 1 . A possibility would be 
that the gauge in these investigations is not the same as in 
the continuum calculations. This will be discussed here. 
Unfortunately, the increase in complexity to fix alterna- 
tive versions of the Landau gauge in the non-perturbative 
domain prevents the access of large lattice volumes at the 
current time in four dimensions. This discrepancy also 
exists in three dimensions [43l . |45| . |47| w hile no discrep- 
ancy exists in two dimensions [45|, S3, [54j . It is therefore 
useful to address the problem first in lower dimensions, 
where large volumes can be accessed more easily. There 
is another advantage when studying the problem first in 
lower dimensions: The fixing of the alternative gauges is 
easier in lower dimensions, as will be shown in section 
IIVI As a result, it will be found that the finite- volume 
effects between both gauges differ. It will also be dis- 
cussed whether they may be the reason for the apparent 
discrepancy. 



The organization of the presentation of this investiga- 
tion is as follows: In section HH general aspects of fixing 
Landau gauges in the non-perturbative domain are dis- 
cussed. Some specific properties of gauge-fixing in con- 
nection with Gribov copies are discussed in section IIIII 
Afterwards, the direct comparison of the propagators in 
the minimal and absolute Landau gauge in two and three 
dimensions is shown in section IIVI Some remarks on the 
four-dimensional case will also be made there. The re- 
sults obtained are discussed in section [V] In particular, 
a possible interpretation will be provided in subsection 
IV Bl Finally, everything is summed up in section IVII 
The algorithm to fix (or, at least, approach) the absolute 
Landau gauge employed here is presented in appendix IB"1 
Technical details of the lattice simulations performed can 
be found in appendix [A"l 



1 See, however, [36, 37] for a possible different interpretation. 



II. LANDAU GAUGES IN THE 
NON-PERTURBATIVE DOMAIN 

In perturbation theory, Landau gauge is defined by the 
local condition 

MS - (i) 

on the gauge field . It belongs to the class of covariant 
gauges, which in general average over the (perturbative) 
gauge orbit with a given weight function [2J. In case of 
the Landau gauge, this weight function is a <5-function. 
Unfortunately, it has been found [H, [6(| that this local 
condition is insufficient to select a unique representative 
of the gauge orbit: There are different elements of the 
gauge orbit, so-called Gribov-Singer copies, which all sat- 
isfy the condition Jl}. An example is, e. g., the instanton 
|6l| . The set of all elements of a gauge orbit satisfying 
condition (frj will be called the residual gauge orbit. It 
has been shown that there is no possibility to specify a 
unique element on the residual gauge orbit with any set 
of additional local conditions [6fj |. 

Further specifications of the gauge will therefore nec- 
essarily be non-local. Any of these should reduce to the 
perturbative Landau gauge |T|) in the perturbative limit. 
This implies that the residual global color symmetry may 
not be broken by the gauge fixing, except when this 
breaking is soft, i. e., vanishes in the perturbative limit. 
As, however, the Kugo-Ojima scenario 0] is based on an 
unbroken global color symmetry, the judicious choice is 
not to break it by the gauge condition. 

It has been argued long ago that the gauge-fixing will 
be complete when the gauge field is further restricted to 
the first Gribov region This region is characterized 
by the semi-positiveness of the Faddeev-Popov operator, 
which in Landau gauge is given by 

-d^d^-gf^A*). (2) 

Unfortunately, this is insufficient [62j. Hence, it is neces- 
sary to further restrict the gauge field. A sufficient pos- 
sibility 2 here is to restrict it to the fundamental modular 
region [62j, which is characterized by the condition that 
the gauge-fixed field configuration absolutely minimizes 

HA) = ^jyj d d xAl(x)A«(x), (3) 

where d is the dimensionality and V is the d-volume. It 
can further be shown that all minima of Q with ((2|) 
positive lie in the first Gribov region ||. Hence, the 
fundamental modular region is enclosed in the Gribov 
region. In fact, both have partly a common boundary 



2 Actually, degenerate minima appear on the boundary of this re- 
gion, which have to be identified 62]. This is resolved implicitly 
by the algorithm employed in the lattice calculations below by 
always choosing exactly one gauge copy for each configuration. 
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Q, on which, e. g., instantons lie [61(. The condition ([3]) 
can be rewritten after double Fourier transformation as 



N„ 



d-l 



D?M ( 4 ) 



where D?* is the gluon propagator and N g the number 
of gluons, 3 for SU(2). Therefore, this condition is equiv- 
alent to minimizing the total trace of the gluon propa- 
gator on the residual gauge orbit. Unfortunately, this is 
an extremalization problem in a high (in the continuum 
infinite) dimensional space. It is therefore NP-hard, sim- 
ilar to spin-glass problems. In particular, the complexity 
increases exponentially with lattice sites. This is one ori- 
gin of the problems in gauge fixing in numerical lattice 
gauge theory calculations. 

From these observations, it is possible to design differ- 
ent types of non-perturbative Landau gauges. One is the 
absolute Landau gauge, defined by finding the absolute 
minimum of ([3]). In the spirit of the perturbative Landau 
gauge, this prescription selects exactly one representative 
on the residual gauge orbit, the one in the fundamental 
modular region. A second possibility is to select an arbi- 
trary representative of the residual gauge orbit, which lies 
inside the first Gribov region. This is the most commonly 
adopted choice in lattice calculations, and is termed the 
minimal Landau gauge, see e. g. [HI]. Of course, other 
possibilities exist. One would be to average over the part 
of the residual gauge orbit which lies in the first Gribov 
region, yielding a Landau-Feynman-like gauge. Another 
would be to average over the complete residual gauge 
orbit. The latter is, however, highly non-trivial, but it 



should be possible to construct 64]. Further ones are 
discussed in sections IIIII and [V] As the origin of field 
configuration space, and thus perturbation theory, is in- 
cluded in all these prescriptions, the correct reduction to 
perturbation theory, in the perturbative limit, is guaran- 
teed. 

Most interesting is whether and, if, to which extent 
results depend on this gauge choice. In particular, mini- 
mal Landau gauge is rather simple to implement in lattice 
simulations [63j , although some caveats will be discussed 
in the next section. On the other hand, absolute Lan- 
dau gauge is an unambiguous definition of the gauge, 
and may therefore be implcmcntable, in particular by 
means of condition also in continuum calculations 
[29| . Hence a comparison of these two gauges is of great 
interest. Such comparisons have been performed repeat- 
edly in the past in four dimensions [iol. 1461. l52l. l57l. l58l . 
However, the results have been rather mixed, and no fi- 
nal result has been obtained so far, although effects seem 
likely, in particular in the case of the ghost propagator. 

Therefore, fixing to the absolute Landau gauge is de- 
sirable. Unfortunately, fixing this gauge becomes expo- 
nentially harder with the number of lattice sites. Hence, 
a natural possibility to assess it is to investigate it in 
lower dimensions, where everything else applies as well, 
in particular the state of the results [1 [E EL El EI [11 
M S M, EEL E3, M, [ML M, HEf ■ This will be done 



here, and it will turn out that this, in fact, is simpler and 
provides currently much more insight. 

It has been argued that for expectation values of oper- 
ators consisting of a product of a finite number of gluon 
field operators the results of the minimal and the abso- 
lute Landau gauge should coincide in the continuum and 
infinite volume limit [ij] 3 . This would apply to the gluon 
propagator, although not to the ghost propagator, which 
depends inversely on the gluon field. There is yet no in- 
dication that this is the case at reachable volumes, and 
also the results presented here do not support this conjec- 
ture. It has been argued that this may be a discretization 
artifact (65j . 

It should be noted that, due to the factor in the 
integral ([4]), the far infrared properties of the gluon prop- 
agator are likely irrelevant for the distinction between the 
minimal and the absolute Landau gauge, as, e. g., an ef- 
fect at 100 MeV can be compensated by a 10 d_1 -times 
smaller effect at 1 GcV momentum: Any difference be- 
tween the gluon propagators in the two gauges in the far 
infrared could be irrelevant to the distinction of the two 
gauges. This applies in particular to the value of the 
gluon propagator at zero momentum, which for a non- 
divergent gluon propagator is not contributing to the in- 
tegral. In addition, due to asymptotic freedom the far 
ultraviolet is again irrelevant for the condition ((4]). 

Furthermore, it should be noted that the statement 
(|U) is true for each configuration individually. The quan- 
tity on the right hand side of ((4]) is the gluon propagator 
on each individual configuration. Nonetheless, the state- 
ment is also true on the average, but the distinction may 
be, and in fact is, only visible at very high statistics due 
to the factor and a difference in ((4]) which is found 
to be of only about 0.3% between the minimal and the 
absolute Landau gauge: For the results presented here, 
the maximization condition cannot be checked explicitly 
for the average gluon propagator due to the statistical 
noise. It is, of course, ensured for each configuration. 

The method to fix the absolute Landau gauge is an 
evolutionary algorithm, which is detailed in appendix [Bl 
Note that similar algorithms have been used previously 
[H, [63] • The minimal Landau gauge will be fixed by an 
adaptive stochastic overrelaxation algorithm [EEI, l63| . 



III. GAUGE-FIXING 

There are a number of peculiarities when it comes to 
gauge fixing on the lattice, which have to be taken into 
account when dealing with Gribov copies. One of them 
is the fact that the complexity of gauge fixing is highly 
dependent on the configuration. A measure of the com- 



3 Because has a unique absolute minimum, by virtue of Q this 
can only occur in the sense of the difference becoming arbitrarily 
small, but an exact equality is never possible. 
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plexity is, e. g., the number of steps an iteration method 
needs to gauge fix a given configuration to the minimal 
Landau gauge. It turns out that this number varies by or- 
ders of magnitudes between different configurations. On 
the other hand, the extreme case of very complex con- 
figurations is rare. There is a relation to the so-called 
exceptional configurations (40l . |55| , and hence the same 
name will be used. Furthermore, it turns out that the 
complexity also depends on the initialization of the iter- 
ative algorithm, i. e., the first guess for the correct gauge 
fixing. It is often possible to reduce the number of nec- 
essary sweeps considerably by choosing other initial con- 
ditions. 

It is therefore tempting to restart the gauge-fixing pro- 
cedure if the situation is exceptional, or even drop the 
corresponding configuration altogether, to improve com- 
putational efficiency. However, in general two different 
restarts correspond to two different Gribov copies, which 
is the observation on which many algorithms to fix to the 
absolute Landau gauge are based on [4(| H2, . This 
may of course lead to a bias when sampling the residual 
gauge orbit, if the complexity of fixing the gauge is, e. g., 
related to the distance to the absolute minimum of © 
along the gauge orbit. That such a bias possibly exists is 
indicated by the fact that a correlation between the value 
of the lowest eigenvalue of the Faddeev-Popov operator 
and the number of gauge fixing sweeps has been observed 
[5^ . The aim in this section is to either support this ob- 
servation by searching for an impact on other observables 
as well, or to show that this is not the case. This will 
be done by performing measurements of observables as a 
function of the number of gauge-fixing sweeps necessary 
to fix to the minimal Landau gauge. Note that any of 
these statements apply only to the system investigated 
here, 40 3 at f3 = 4.2, and may be different for different 
systems, although empirically these results seem to be 
typical. 

To investigate this, configurations have been obtained 
from the Wilson action according to [5 a ] and th en g auge- 
fixed to (minimal) Landau gauge according to [55] , The 
details of the configurations are given in table U in ap- 
pendix [X] The gauge-fixing procedure was modified by 
imposing a limiting number of sweeps n s . If during 
gauge-fixing more than this limit would be necessary, 
then the gauge-fixing process is restarted with another 
initialization. The latter was always a randomly cho- 
sen gauge transformation. If it was not possible to fix 
the gauge after twenty restarts, the configuration was 
dropped, and a new one was generated instead. 

The results have been compared to the case without 
any limitation, but otherwise the same gauge-fixing pro- 
cedure. In this case, the maximum number of gauge- 
fixing sweeps observed was n s = 28740, while the av- 
erage number was 1528(55), and the minimum number 4 
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FIG. 1: Top: Fraction of restarts of the maximum possible 
number of restarts during gauge- fixing as a function of n s , 
not counting dropped configurations. Bottom: Fraction of 
dropped configurations compared to the total number of con- 
figurations as a function of n 3 . Values greater than 100% are 
possible, if more than every second configuration has been 
dropped. 



4 This number was due to an optimization of the gauge-fixing pro- 



cedure, and is actually only an upper limit for the minimum 
number. 
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401. To quantify the amount of exceptional configura- 
tions with respect to the limit n s the number of neces- 
sary restarts and of dropped configurations 5 are shown 
in figure [TJ It is visible that very extreme configura- 
tions are very rare. Within an order of magnitude of the 
maximum of gauge-fixing sweeps, it was always possible 
to find another Gribov copy which could be fixed suffi- 
ciently fast. Therefore, severe restrictions on the number 
of gauge-fixing sweeps are necessary for configurations to 
be dropped. Still, starting from a certain limit on, a very 
strong rise in both the number of restarts and the number 
of dropped configurations is observed. 

It turns out that the exceptional configurations have 
no statistically significant impact on observables as long 
as the number of restarts is sufficiently small. As an 
example of the impact on measured quantities, the gluon 
propagator 

Df v (x-y)=<A«(x)A» u (x)> 

and the ghost propagator 

D^ix -y) = -< (d,(d,5 ab - gf^AUx)))- 1 > 

will be used. In both cases, the color (and, for the gluon 
propagator, Lorentz) traced 6 and Fourier-transformed 
propagators will always be shown, D(p) and Dg{p), re- 
spectively 7 . 

The propagators are shown in figure[2j It is visible that 
there is essentially no change for the gluon propagator. 
In case of the ghost propagator the few lowest momentum 
points are affected, which would be in accordance with 
the observation of an influence on the lowest eigenvalue 
of the Faddeev-Popov operator [55[ . This latter influence 
is also observed here, manifesting itself by a systematic 
increase of the lowest eigenvalue of the Faddeev-Popov 
operator by about 30% 8 . However, the lowest two mo- 
mentum points of the ghost propagator are subject to 
finite volume effects [53l [56j . and the consequences are 
thus less relevant. Thus, at the respective level of biasing 
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FIG. 2: The gluon propagator (top panel) and the ghost dress- 
ing function (bottom panel) as a function of n s for n a = oo 
(full circles), n s = 2212 (open circles) and n s = 791 (stars). 



5 Indications have been found that it is very likely that if one con- 
figuration is dropped the next one(s) will be dropped also. This 
would indicate that the correlation time increases for configura- 
tions which are harder to gauge-fix, as the number of sweeps 
between two gauge-fixed measurements is sufficient for many 
other observables to decorrelate. To make this a more reliable 
statement more investigation is needed, however |68l| . To com- 
pensate this effect, multiple independent runs have always been 
performed throughout all results in this paper. 

6 None of the modifications of the gauge-fixing process performed 
here or in section llVl changes the fact that both propagators are 
color-diagonal lo4l. lo r l . 

7 A more precise definition can be found in [55], along with details 
on the determination of these quantities, including error deter- 
mination. 

8 The presented deformation of the sampling procedure of the 
residual gauge orbit would therefore influence the limits on the 
infinite volume ghost propagator constructed in |47|| , if it would 
be determined with the limits presented here. 



the selection of Gribov copies has almost no consequence 
within statistical errors for the propagators. As a further 
observable the gluon propagator at zero momentum will 
be investigated. This value is particularly important, as 
will be discussed in the next section. Its dependence on 
n s is shown in figure [31 Again, within statistical errors, 
there is no consequence of limiting the number of gauge- 
fixing sweeps. 

Therefore, at least for a modest restart policy and at 
the given level of statistical errors, for the system inves- 
tigated here, the consequences of the bias observed in 
|55| have no relevant impact, except for the properties 
of finite volume effects. However, as will be seen in the 
next section, a bias in selecting Gribov copies can have 
consequences. Therefore, when desiring to use a limita- 
tion to improve the usage of computational resources, it 
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FIG. 3: The gluon propagator at zero momentum as a func- 
tion of the gauge fixing sweeps limit n„ . 



would be necessary to make sure in advance whether for 
a given system and a given observable this introduces a 
bias. Hence, no such selection criteria will be used here- 
after. 

The reason for this effect may be that Gribov copies 
which are more simple to fix to the Landau gauge are 
in some sense more 'perturbative', which is by defini- 
tion almost in Landau gauge. This would also explain 
the increase of the value of the lowest eigenvalue of the 
Faddeev-Popov operator. This is also supported by an- 
other observation: When the thermalization process is 
initialized with a cold, i. e. vacuum configuration, and 
not enough thermalization sweeps are performed to fully 
thermalize the configuration, it is usually much simpler 
to gauge fix than a fully thermalized configuration. 

However, as there is an effect on correlation functions 
by performing this modification of the minimal Landau 
gauge it implies that 'minimal Landau gauge' is not by it- 
self a sufficiently precise description of the gauge. Implic- 
itly assumed is that dropping configurations is only a con- 
venient numerical simplification, and is not altering the 
sampling of configurations. This is supported by the fact 
that the alterations observed have a smooth onset even 
when none or only a statistically insignificant number of 
configurations are dropped. A more precise description is 
hence needed to qualify what is meant by 'choosing ran- 
domly one representative inside the first Gribov region'. 
This has to be understood further to give a precise mean- 
ing to the concept of the minimal Landau gauge. Here, as 
an operational definition of the gauge the description of 
the gauge-fixing algorithm (with no limit on the number 
of sweeps) [55[ will be used. 

This point will be discussed further in section fVl 



IV. MINIMAL VS. ABSOLUTE LANDAU 
GAUGE 

In the following the propagators, as examples of gauge- 
dependent Green's functions, in minimal and absolute 
Landau gauge will be compared. Aside from the gen- 
eral question to which extent gauge-dependent correla- 
tion functions actually differ in the two gauges, there is 
another reason to pose this question. Predictions exist 
from functional calculations on the infrared behavior of 
these correlation functions 0, H, HH, E3 j m fact, for any 
correlation function [IH, . These predict a power-law 
behavior of the correlation-functions in the far infrared, 
which is, if it exists, unique [16l |. In particular, the ghost 
propagator would be more singular than that of a mass- 
less particle, and the gluon propagator would be infrared 
vanishing if the ghost-gluon vertex were infrared con- 
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The latter is, inside this scenario, 
and supported by lattice calculations 



5Ct l54[ . This applies to two, three and four dimcn- 
It has been discussed, whether this could be the 
form of the propagators in the absolute Landau gauge 

However, as another possibility it has been suggested 
that the gluon propagator should be infrared finite and 
the ghost propagator behave like the one of a massless 
particle. This has been pro posed in three and four di- 
mensions [n m, m m, in m, i3, n m, m m, 113 . 

In particular, this has been obtained in the context of a 
minimal Landau gauge [36|, H3 • 

The functional methods generally suffer from necessary 
assumptions and approximations to make them solvable. 
Although this necessity has been greatly reduced over 
time OIHjISIlll) the question which is the infrared be- 
havior of the correlation functions still induces uncertain- 
ties. Therefore, it is natural to compare these results to 
lattice calculations, which have no need for assumptions. 
On the other hand, they suffer from artifacts due to finite 
volumes, discretization and numerical requirements. In 
particular, finite volumes have been a serious issue, espe- 
cially in four dimensions [53]. Recently, however, compu- 
tational resources have become available, which reduced 
this problem considerable, at least for Yang-Mills theory 
with sufficiently coarse discretization [H, |3, EE 53 ■ 

A serious obstacle is still the necessity to fix to the 
absolute Landau gauge. This is a minimization problem 
in a high-dimensional space, therefore NP-complete, and 
thus exponentially hard. However, arguments have been 
made that at least for the gluon propagator in the limit 
of infinite volume in the continuum the results should be 
identical in the minimal and the absolute Landau gauge 
U Therefore, the numerically much simpler, since poly- 
nominally hard, minimal Landau gauge has been used to 
obtain results for the propagators on the lattice. These 
are in good agreement with an infrared finite gluon prop- 
agator and a massless ghost propagator in four and three 
dimensions [H [44], Ea, S3, but an infrared vanishing 
gluon propagator and infrared enhanced ghost propaga- 
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tor in two dimensions [45j, |47J, |54| . 

Still, the claim that both gauges coincide on the level 
of propagators has to be verified. Various attempts have 
been made, all on (necessarily) small volumes in four di- 
mensions IH, EH, [57], [H, H^]. Some discrepancies 
have been found for the ghost propagator between both 
gauges, which may or may not decrease with volume. 
However, as the problem increases with the dimensional- 
ity of the space, and thus lattice sites, it can be expected 
that it is much simpler in lower dimensions. Furthermore, 
it has been found that in two dimensions the results are 
in favor of a power-law solution [54j . It is therefore sug- 
gestive that also in higher dimensions in fact the power 
law solution would persist in the absolute Landau gauge. 
This is supported by the fact that in three dimensions, 
where the problem should be more severe than in two 
but less than in four, at least partly a power-law solution 
is observed [5l|. and a discrepancy only appears at very 
large volumes [45j, [47J . It is therefore a natural investiga- 
tion to check the situation in two and three dimensions. 
This will be done here. 

To this end, an evolutionary algorithm [69j will be 
used, which is described in detail in appendix [Bj Still, 
this algorithm is not able to guarantee to fix to the abso- 
lute Landau gauge, as it cannot guarantee to find the ab- 
solute minimum of the gauge-fixing functional ([5]). This 
is a general problem gfl, SI, [H [53, [H, HE Is3 for this 
type of minimization problem, which belongs to the same 
class as finding the ground state energy of spin-glass mod- 
els. Only in the asymptotic limit of infinite computing 
time success can be guaranteed. 

This has the unpleasant consequence that, even if the 
algorithm actually succeeded in finding the absolute min- 
imum, it would be impossible to tell so. The only thing 
which can be hoped for is that the found gauge-fixed con- 
figuration is 'closer' to the absolute Landau gauge. Even 
a notion of 'closer' is impossible to define rigorously, as it 
is even unknown whether any type of gauge-fixing of the 
residual gauge orbit is smoothly connected to the abso- 
lute Landau gauge or not. That this is the case will be 
assumed henceforth. 

Therefore, as a pragmatic ansatz, the fulfillment of the 
condition ((4]) , finding the gluon propagator with the min- 
imum total trace, will be taken as a criterion to compare 
the closeness to the absolute Landau gauge of two gauge- 
fixed configurations. Hence, the gauge-fixed configura- 
tion with the least value of J-(A) found will be taken 
to be in the absolute Landau gauge. In addition, as a 
measure of the presence of Gribov copies the difference 



1 A as a function of volume I 



A = T(A worst ) " F(A hest ) 



(5) 



will be used, where 'best' and 'worst' refer to the gauge- 
fixed field for a given configuration having the lowest and 
highest value of !F{ A) , respectively. 

It turns out that the function ([5]) has a generic shape, 
with potentially three, rather sharp, peaks. The lowest 
has values of A of the order of the numerical precision, 
and is therefore clearly due to numerical noise. The sec- 
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FIG. 4: The presence of Gribov copies as measured by A as 
a function of inverse volume in two dimensions. Full circles 
correspond to (3 = 10, while open circles correspond to (3 = 
38.7. 



ond peak clusters at very small values, typical of around 
the order of 10 -15 , likely due to lattice Gribov copies 
f70j |. The third peak is around 1CP 3 and appears to be 
related to continuum Gribov copies. With increasing vol- 
ume, the population moves from the lowest to the highest 
peak. At sufficiently large volumes, only the peak at the 
highest value of A is populated. This could indicate that 
then at least one continuum Gribov copy is present in 
each configuration, in addition to any lattice copies. 

Technical details on the lattice simulations performed 
can be found in appendix [A"l 



A. Two Dimensions 

The case of two dimensions is particularly interesting, 
as the results in lattice calculations and in continuum 
functional calculations agree @, [54[. Furthermore, in 
this case the finite-volume effects for the gluon propa- 
gator behave in exactly the same way as in functional 
calculations, while they are apparently different for the 
ghost propagator 0, Ell . It is therefore highly interest- 
ing, whether a difference in terms of correlation functions 
between the minimal and absolute Landau gauge exists, 
and, in particular, whether this will harm the agreement 
between the continuum and the lattice results. 

The dependence of A on the (inverse) volume is shown 
in figure 2J Even at quite large volumes, A is still rather 
small, implying that a considerable amount of the config- 
urations do not yet turn up in the peak with the highest 
value of A. This would indicate that in many cases the 
difference between the best and the worst Gribov copy 
is still negligible in terms of the gluon fields themselves. 
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FIG. 5: The propagators from a 136 2 ((25 fm) 2 ) lattice at p = 
10. The top panel shows the gluon propagator and the bottom 
panel the ghost dressing function. Full circles correspond to 
the results in the minimal Landau gauge and open circles to 
the absolute Landau gauge. 



On the other hand, the exponential rise with volume is 
still clearly visible. It cannot, however, be anticipated at 
which level the function finally levels, and thus, whether 
Gribov copies could have a sizeable impact on the gluon 
fields at sufficiently large volumes. This remains to be 
seen. Furthermore, the discretization has only a small 
effect. In addition, this effect decreases with increasing 
volume. 

However, a statement on the average difference be- 
tween the field variables is not sufficient for a statement 
on their correlation functions. Therefore, these are stud- 
ied separately. For a very large volume, the difference 
between the propagators in minimal Landau gauge and 
in absolute Landau gauge are shown in figure [5J The im- 
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FIG. 6: The gluon propagator at zero momentum as a func- 
tion of inverse volume. Full and open symbols correspond to 
the minimal and absolute Landau gauge, respectively. Circles 
are from calculations at /3 = 10 and squares from calculations 
at p = 38.7. 



pact is quantitatively rather small for the propagators in 
the two different gauges. A statistically significant effect 
is only seen for the ghost propagator at momenta below 
250 MeV. There, the propagator in the absolute Landau 
gauge is systematically smaller than in minimal Landau 
gauge. However, the effect is strongest for the two lowest 
momentum points, which are anyhow unreliable due to 
finite volume effects. There is also a small effect on the 
border of statistical significance for the gluon propaga- 
tor. The gluon propagator in absolute Landau gauge is 
smaller in the infrared than the one in minimal Landau 
gauge. This would be in agreement with the condition 
((3]). Nonetheless, the strongest effects are in the domain 
where finite volume artifacts dominate. Therefore, the 
effect is almost negligible. Still, the volume dependence 
of the gluon propagator at zero momentum is a quantity 
for which predictions exist 15 611 . and which in the minimal 
Landau gauge are fulfilled [5J]. In particular, it behaves 
as a power-law and vanishes in the infinite- volume limit. 
This is not altered in the absolute Landau gauge, as seen 
in figure [5] As can be seen, there is essentially no differ- 
ence in the volume dependence. 

A final possibility for comparison are differential char- 
acteristics of the propagators. Both propagators are pre- 
dicted to behave as power-laws in the far infrared region 
0, [ll| . This can be parameterized by 

D{p) « P «A YM Azip 2 )- 1 -' 
D G (p) « p «Avm Aaip 2 )- 1 -*, 

with constant prefactors Az and Aq. In particular, the 
effective infrared exponents t and function of vol- 

ume are interesting, as they are the key quantities in 
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FIG. 7: The gluon exponent t (top panel) and the ghost ex- 
ponent k (bottom panel) as a function of volume. Full and 
open symbols correspond to the minimal and absolute Landau 
gauge, respectively. Circles are from calculations at /3 = 10 
and squares from calculations at [3 — 38.7. The lines indicate 
the infinite-volume predictions from functional calculations 
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comparison to functional calculations. This volume de- 
pendence of the exponents, determined in the same man- 
ner as in [54| , is shown in figure [7] The gluon exponent t 
does not differ in both gauges, and in both cases it tends 
towards the predicted value in the infinite volume limit. 
In case of the ghost exponent, there seems to be a sys- 
tematic difference, both in central value and in volume 
dependence between both gauges. However, this fact is 



Minimal Landau gauge 
Absolute Landau gauge 




FIG. 8: The quantity © from a 136 2 ((25 fm) 2 ) lattice at 
13 — 10. Full circles correspond to the results in the minimal 
Landau gauge and open circles to the absolute Landau gauge. 
Note that the two lowest momentum points are affected by 
finite volume artifacts [54j |. 



not statistically significant 9 . In addition, in both gauges 
the exponent is compatible to the predicted value, within 
statistical accuracy. The variations of the central values 
for the two largest volumes indicate that the results at 
the respective level of statistics become already unreli- 
able, a consequence of the fit procedure [54j . 

A more sensitive test is the relation predicted between 
both exponents [|, 1, [H 13] 



t = -2k 



4-d 



As a consequence, the quantity 

a{p) = D(p)D G (p) 2 p 2+d 



(6) 



(7) 



should be a constant in the far infrared. This quantity 
is proportional to a possible definition of the coupling 
constant [ll| and can be regarded as an effective coupling 
constant. 

The relation © holds in both gauges, as can be seen 
from the quantity ([7]) shown in figure [8l However, the 
constant of proportionality shows a difference between 
both gauges. This implies that the difference observed 
for the ghost propagator in both gauges in figure [5] is 
likely rather due to a change in the prefactor Aq from 
gauge to gauge than in the value of the exponent k. 

Hence, it seems that there is at best a quantitative 
difference between the propagators in both gauges, but 



9 With the statistics employed in 
have been visible at all. 



this difference would not 
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FIG. 9: The presence of Gribov copies as measured by A 
as a function of inverse volume in three dimensions. Full 
circles correspond to (5 = 4.24, while open circles correspond 
to p = 7.09. 



FIG. 10: The presence of Gribov copies as measured by A as 
a function of inverse volume in three dimensions, averaged on 
a linear scale. Full circles correspond to f3 = 4.24, while open 
circles correspond to f3 = 7.09. 



not a qualitative one. Even this quantitative difference 
seems to be rather small. In particular, the propagators 
in both gauges agree with the predictions from functional 
calculations in the continuum. However, after investigat- 
ing the case in three dimensions in the next section, it 
will be necessary to reconsider this conclusion. This will 
be done in section IV Bl 



B. Three Dimensions 

In three dimensions, the situation is quite different 
than in two dimensions. The predictions from functional 
calculations are qualitatively identical [f| [ill, US H3| • Re- 
sults on small and medium-sized lattices seem to confirm 
these predictions [HTI . l55l [7l| . However, results on the 
largest lattices suggest rather an infrared finite gluon 
propa gato r and an infrared tree-level-like ghost propa- 
gator [45l H3| , as obtained in functional calculations in 
[37| . Still, at least in an intermediate momentum win- 
dow these results exhibit a power-law-like behavior. The 
latter observation is not a finite volume effect. 

Starting again with the measure for the presence of 
Gribov copies A, the result is already drastically different 
from the two-dimensional case. This is shown in figure 
[9l As in the two-dimensional case, an exponential rise 
with volume is visible. However, it saturates at a value 
of about 10~ 3 , corresponding to the highest-valued peak 
in the distribution of A. In fact, at the largest volumes, 
A is confined to a rather small region around this peak, 
implying that in all configurations there are significant 
differences in terms of ^F(A) between the best and worst 
Gribov copy. If these large differences in fact correspond 



to continuum Gribov copies, this implies that, starting 
from volumes of the order of about (3 fm) 3 , Gribov copies 
play an important role. 

A little detail is that A seems to decrease once more 
when the volume becomes larger. This is more drasti- 
cally seen when averaging A on a linear scale, as shown 
in figure [TUJ Here, it is visible that A starts to decay 
significantly with increasing volume. On the one hand, 
this can indicate that the local minima become more and 
more degenerate with the absolute minimum, and is thus 
a statement on the nature of Gribov copies 10 . It could 
also imply that the gauge-fixing algorithm becomes un- 
able to find a good approximation of the absolute min- 
imum. This would be possible if the search space cov- 
ered is not sufficiently enlarged to compensate for the ex- 
ponential increase of the number of Gribov copies with 
volume. This possibility is discussed in more detail in 
appendix [B] and is the most likely one. It could, of 
course, also be a mixture of both effects. Anyway, it 
may potentially be affecting the reliability of the results 
at sufficiently large volumes. Eventually, if the compu- 
tational resources are not increased exponentially with 
volume as well, the gauge-fixing algorithm will fail. As 
a consequence, the result will become indistinguishable 
from the ones in minimal Landau gauge. This may be 
the onset of it. A third possibility is that the Zwanziger 
conjecture may be correct in the end: A decrease in A 
implies that the integrated strength of the gluon propa- 



This happens, e. g., in 2d U(l) gauge theory |72j . I am grateful 
to L. von Smekal for pointing out this possibility. 
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gator becomes more and more similar between minimal 
and absolute Landau gauge once more, leading possibly 
eventually to coinciding gluon propagators. However, it 
will be shown below that this is not equivalent to a point- 
wise coincidence of the gluon propagators yet, although 
some degradation is observed for the results in the largest 
volume. 

Monitoring the gluon propagator is another possibil- 
ity to measure the success of the gauge-fixing procedure. 
The result for the gluon propagator in both gauges is 
shown in one large volume in figure llll Alongside also 
the volume dependence of the propagator in the infrared 
and the behavior of it at zero momentum are shown. 
First of all, the statistically significant difference is quan- 
titatively rather small, and confined to momenta below 
roughly half a GeV. The latter is expected, as the dif- 
ference between both gauges is purely non-perturbative, 
and should therefore decrease as a power of momentum 
towards large momenta, leaving only the identical per- 
turbative part of the propagators. 

Interestingly, the statistically significant deviations 
lead to a gluon propagator which is smaller at low mo- 
menta than the one in minimal Landau gauge, in explicit 
agreement with condition Any possibly existing de- 
viations at larger momenta are not visible due to the 
statistical fluctuations 11 . 

Hence, the gluon propagator in absolute Landau gauge 
is more strongly infrared suppressed than in minimal 
Landau gauge. In particular, the gluon propagator at 
zero momentum decays faster in absolute Landau gauge 
than in minimal Landau gauge, and continues longer in a 
power-law decay towards the infinite- volume limit, while 
in minimal Landau gauge the gluon propagator seems to 
exhibit a mass-like behavior, as has been observed previ- 
ously [H|,[45|]. Only for the largest volume the value of the 
propagator at zero momentum seems to rise once more. 
However, as there is no similar effect for non-zero momen- 
tum, and the statistical error is comparatively large, this 
may, or may not, be a statistical fluctuation. Otherwise, 
it may be an effect due to the one observed in figure [TU] 
that A decreases for large volumes. Beyond this, at least 
for the volumes investigated here, this decrease does not 
(yet) have an effect. In particular, excluding the largest 
volume, the value of D(0) in minimal Landau gauge can- 
not be fitted with a pure power-law, while this is possible 
in the absolute Landau gauge. This is shown in figure 
llll In general, the gluon propagator in absolute Landau 
gauge therefore could be vanishing at zero momentum in 
an infinite volume. In minimal Landau gauge, however, 
it seems to develop a mass, and becomes infrared con- 
stant. Furthermore, both behave similarly but not the 



11 Note that the results for both gauges are obtained from two 
independent sets of configurations to reduce systematic errors. 
Hence a direct configuration-by-configuration comparison is not 
possible. Furthermore, fulfillment of the condition (|4jl may be 
implemented at different momenta for different configurations. 



same as a function of volume. However, larger volumes 
and better statistics will be necessary to substantiate this 
observation. 

Corresponding results for the ghost propagator are 
shown in figure 1121 The results are highly interesting. 
Comparing within a small volume the propagator in both 
gauges, it is found that the one in absolute Landau gauge 
is significantly smaller up to momenta as large as nearly 1 
GeV than the one in minimal Landau gauge. In a larger 
volume, however, both look much more similar again. 
The reason for this behavior can be seen from the volume- 
dependence of the propagator in both gauges. In mini- 
mal Landau gauge, the propagator becomes less and less 
infrared enhanced with increasing volume, until it even- 
tually becomes essentially the one of a massless particle 
[43l \4l\ . This is, however, not the case in the absolute 
Landau gauge. While it also becomes less divergent for 
small volumes, the situation changes when Gribov copies 
become relevant. Then, the propagator stays at least as 
divergent as it was on smaller volumes, or may possibly 
become even more divergent again. 

These results are reflected also in the spectrum of 
the Faddeev-Popov operator, in particular on its lowest 
eigenvalue. However, the statistical uncertainty in this 
case is too large for a definite conclusion. 

To complete the collection of results on the propaga- 
tors, two more properties are investigated. One is again 
the differential properties in the form of the effective ex- 
ponents t and k, which are shown in figure [T3l The expo- 
nent for the gluon propagator is mildly affected at large 
volumes, becoming more negative in the absolute Landau 
gauge. Still, the volumes employed here are such that 
also in minimal Landau gauge the exponent is smaller 
than —1, but larger than the one in absolute Landau 
gauge. Both are thus for these volumes in accordance 
with an infrared vanishing gluon propagator. Of course, 
this is no longer the case on larger volumes in case of 
the minimal Landau gauge [43| . but it is currently un- 
known what will be the case for the absolute Landau 
gauge. In case of the ghost exponent, the difference be- 
tween both gauges is more pronounced, but statistical 
uncertainties prevent a definite conclusion at large vol- 
umes. In particular, as can already be deduced from the 
ghost dressing function in figure [121 the ghost exponent 
is initially larger in minimal Landau gauge than in ab- 
solute Landau gauge. The exponent in minimal Landau 
gauge, however, decreases with volume [43l |47| while the 
one in the absolute Landau gauge is essentially volume 
independent within statistical accuracy. This once more 
demonstrates that the results in absolute Landau gauge 
exhibit different finite volume effects. 

Finally, the quantity a, defined in equation ([7]), should 
become constant for all dimensions in the far infrared, 
and thus also in three dimensions. This is the case in 
two dimensions in both gauges, as shown above. How- 
ever, it seems to go to zero in three and higher dimensions 
[43| . This is to be expected in any finite volume [56| , but 
the pattern of finite volume artifacts observed in mini- 
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FIG. 11: The gluon propagator in three dimensions. The top left panel shows the comparison between minimal (full circles) 
and absolute (open circles) Landau gauge on a 56 3 lattice at f3 = 4.24, corresponding to V = (9.5 fm) 3 . The bottom panels 
show the volume dependence for the gluon propagator in minimal (left panel) and absolute (right panel) Landau gauge. Circles, 
squares, and triangles correspond to V = (11 fm) 3 , (5.4 fm) 3 , and (2.7 fm) 3 , respectively, all calculated at j3 = 4.24. The top 
right panel shows the gluon propagator at zero momentum, together with fits for the coarser discretization. Full symbols are 
from the minimal Landau gauge, while open symbols are from the absolute Landau gauge. Circles and squares correspond to 
calculations at (3 = 4.24 and /? = 7.09. 



mal Landau gauge is not in agreement with a constant 
value in the infinite volume limit [ill l47j ] . The result in 
absolute Landau gauge, and its volume dependence, is 
shown in figure [141 The difference between both gauges 
is not qualitative, in particular, in both cases a seems 
to go to zero. However, there is an interesting volume 
dependence observed in the absolute Landau gauge. It 
is seen that with increasing volume the maximum tends 
to move towards smaller momenta, albeit rather slowly. 
The consequence of this is particularly visible in case of 
the second- lowest momentum point: The same value of 



a is attained at smaller and smaller momenta on the de- 
creasing slope. This would be in accordance with the 
type of finite volume effects predicted for an infrared fi- 
nite coupling constant in the infinite- volume limit [56| . 



C. A note on four dimensions 

Take now for a moment for granted that in absolute 
Landau gauge the infrared behavior of the Green's func- 
tions would be given by power laws. This would imply 
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that the effective gluon exponent would become smaller 
with increasing dimension, while the one for the ghost be- 
comes larger [5|. Given the volumes from which onward 
a maximum in the gluon propagator can be observed in 
two and three dimensions, it can be immediately deduced 
that lattice volumes of the order of 50 4 at (3 — 2.2 would 
be necessary to obtain such a maximum unambiguously 
in four dimensions. Furthermore, from the results for 
A obtained on small lattices, which is shown in figure 
IT51 and should be compared to the figures [4] and 02 it is 
clear that the number of Gribov copies in four dimen- 
sions at a given linear lattice extension in lattice units 
significantly exceeds the same number in lower dimen- 



sions. As assured in the beginning the issue of Gribov 
copies becomes more severe with an increasing number 
of dimensions. 

From this already a rough estimate of the necessary 
computing time can be deduced, which will be needed 
to obtain the propagators in absolute Landau gauge for 
the indicated volume. This exceeds the amount of com- 
puting time invested here for the the calculations in two 
and three dimensions by more than an order of magni- 
tude. Aside from this, the necessary amount of storage 
space at intermediate steps of the calculations to hold all 
the gauge transformations needed in parallel for the algo- 
rithm presented in appendix IBl is also rather prohibitive. 
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FIG. 13: The gluon exponent t (top panel) and the ghost 
exponent n (bottom panel) as a function of volume in three 
dimensions. Full and open symbols correspond to the minimal 
and absolute Landau gauge, respectively. Circles are from 
calculations at /3 = 4.24 and squares from calculations at 
f3 = 7.09. 



Hence, this has to be investigated in detail further, when 
sufficient resources will be available in the near future. 



V. DISCUSSION 

A. Non-perturbative definitions of the Landau 
gauge 

Within the first Gribov region many Gribov copies ex- 
ist. All of the gauge-fixing techniques used in this work 
select for each configuration one particular of these min- 
ima. However, which of these minima is selected depends 



FIG. 14: In the top panel the quantity a, defined in J?} is 
shown in minimal Landau gauge (full circles) and in absolute 
Landau gauge (open circles). In the lower panel the volume 
dependence of the effective coupling constant in absolute Lan- 
dau gauge is shown. Circles are from a 64 3 lattice, squares 
from a 56 3 lattice, and triangles from a 48 3 lattice, all calcu- 
lated at fi = 4.24. 



on the algorithm. Two different versions which use an 
implicit probabilistic choice are the adaptive stochastic 
overrelaxation procedure, which has been used as defin- 
ing the minimal Landau gauge throughout, with and 
without a limit on the number of gauge-fixing sweeps. 
The probabilistic aspect comes into play as the Gribov 
copy chosen depends on the initial (random) guess for 
the correct gauge transformation to Landau gauge. In 
this sense, both lead to different average values of the 
gauge-fixing functional J- (A) ([3]), defined by their differ- 
ent weighting functions. Although the differences on the 
level of Green's functions are small, there exist differences 
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FIG. 15: The presence of Gribov copies as measured by A as 
a function of inverse volume in four dimensions. Full circles 
correspond to f3 = 2.2, while open circles correspond to J3 = 
2.45. 



for both possibilities, as shown explicitly in section ITTT1 

An alternative definition is by the selection of the low- 
est minimum (found) . This is the method used for the ab- 
solute Landau gauge. Again, this corresponds to a differ- 
ent value of !F(A), and also to different Green's functions, 
as has been shown in section IIVI It would, of course, 
be possible to select instead of the lowest minimum, the 
second-lowest, the third- lowest etc., which gives a one- 
parameter-family of Landau gauges 12 . Of course, all of 
them, by construction, lead to a different value of ^(A), 
and therefore, by virtue of to differing gluon propa- 
gators. 

This, so far, results from the possibilities how to de- 
fine the Landau gauge, and are consequences of the non- 
perturbative ambiguities of the Landau gauge condition 
©■ 



B. A possible interpretation 

It is now possible to interpret these concepts and re- 
sults using an existing framework [29] . As the data is not 
yet yielding unambiguously one result, given the limited 
resources available, this is, of course, just a possibility. 

It is conceivable that gauges based on selecting a par- 
ticular minimum out of the list of ordered minima will 
yield Green's functions with values between the results 



12 On a finite lattice, there will be no degeneracies. In the infinite- 
volume and continuum limit, this may be different, and a pre- 
scription how to deal with these degeneracies would be required 
as well. 



in the minimal Landau gauge and the absolute Landau 
gauge, as long as the respective values of F(A) are be- 
tween those of the absolute and the minimal Landau 
gauge. In case of values larger than the one of minimal 
Landau gauge, it is not a-priori clear what might occur. 
Although this conjecture is likely, it requires explicit fu- 
ture proof. 

In particular, it would be interesting to track whether 
the value of -D(O) or p 2 Da{p)\ p =o could be mapped to 
such a definition of the non-perturbative Landau gauge. 
If this would be the case, it is tempting to assume that 
this family of Landau gauges could be in one-to-one corre- 
spondence to the one-parameter family of solutions found 
for the corresponding functional equations 29] . Further- 
more, this would give a direct explanation why the cor- 
responding effective action in minimal Landau gauge is 
not having a unique minimum (to one-loop order) as a 
(implicit) function of D(0) [33], at least in three dimen- 
sions: The value of D(0) (or better p 2 D G (p)\ p=0 
would then have no physical significance, but only acts 
as a gauge parameter. This would be, of course, highly 
welcome, as D(0) cannot be an independent parameter 
of the theory, which is completely specified once the cou- 
pling, or the scale, has been fixed. On the lattice, of 
course, this family is finite, as only a finite number of 
Gribov copies exists (or at least are found), while this is 
a continuous one-parameter family in the continuum. 

A possibility would then be that the absolute Landau 
gauge in fact corresponds to the endpoint of this one- 
parameter family, i. e., to the value (p 2 Dc(p))~ 1 \ P =o = 0. 
Whether in the continuum the scaling solution possibly 
fulfills the condition Q for the absolute Landau gauge, 
and to which extent this affects solutions with a finite 
value for (p 2 Dg(p))~ 1 \ p =o is discussed in [29j. However, 
these are currently truncation-dependent statements [291 ] . 

Also the results obtained in Landau gauge here ex- 
hibit a dependence on volume, which would be compat- 
ible with the one expected for a power-law solution in 
the infinite- volume and continuum limit [561 ] . although 
the results are not unambiguously in favor of this, they 
are neither in strict disagreement. In particular, all 
changes which occur when moving from the minimal Lan- 
dau gauge to the absolute Landau gauge are exactly as 
expected from the scenario sketched here. 

The only, at first sight, disagreeing result is the ghost 
propagator in absolute Landau gauge being smaller than 
in minimal Landau gauge. The reason for this behav- 
ior of the ghost propagator, shown in figure Q2] in three 
dimensions, may be connected to the volume evolution 
of the first Gribov region and the fundamental modu- 
lar region. Due to volume effects, in both cases config- 
urations at the boundary dominate the average. How- 
ever, at any finite volume both do not have a common 
boundary ||. Given that the first Gribov region is 
bounded by configurations with the smallest eigenvalues 
of the Faddeev-Popov operator, it appears that the ghost 
propagator has to diverge less in absolute Landau gauge 
than in minimal Landau gauge [4?], l59j . But not only the 
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eigenvalues enter the determination of the ghost propa- 
gator, but also the eigenstates [73[. Although occasion- 
ally assumed, the Faddeev-Popov operator does not have 
in the infinite-volume and continuum limit only normal- 
izable eigenstates. This is due to its equivalence to a 
quantum mechanical Hamilton operator, which can have 
scattering states as eigenstates. In particular, explicit ex- 
amples exist for cases where the lowest (zero) eigenvalue 
has a non-normalizable eigenstate [6l|. This is also an 
important fact in the Gribov-Zwanziger framework 5, 6]. 
On the other hand, in any finite volume, the eigenstates 
can be normalized, and taking the (possible) divergence 
of this normalization in the infinite- volume limit into ac- 
count could be crucial. The combined effect of the pre- 
cise shape of the low- eigenvalue spectrum and a diverg- 
ing normalization of the eigenstates could be sufficient to 
provide a more infrared divergent ghost propagator in the 
infinite-volume and continuum limit in absolute Landau 
gauge than in minimal Landau gauge. Identifying how 
this could actually be realized is, however, beyond the 
scope of this work. A final proof is, of course, only that 
in the absolute Landau gauge at sufficiently large volume 
the ghost propagator would be more singular than in the 
minimal Landau gauge. The volume-dependence of the 
propagator in both gauges found here is as expected if 
this is the case. 

A consequence of this scenario is that it should be ex- 
pected that also in two dimensions, for sufficiently large 
volumes and number of Gribov copies, an infrared fi- 
nite gluon propagator is obtained in the minimal Landau 
gauge. Given the results shown in figure HI this must be 
at volumes much larger than the ones currently investi- 
gated. This will be checked in the future. Correspond- 
ing observations have already been made analytically in 
Coulomb gauge in two dimensions [74j . 

Finally, the question of discretizations artifacts has so 
far not been sufficiently addressed, although indications 
exist that these might still be relevant [63. This will 
also be investigated in the future. In other gauges, e. g. 
linear covariant ones, they are of considerable importance 

MM- 

In addition, from the lattice perspective it is not fully 
solved how this relates to the Kugo-Ojima scenario, as 
this requires the construction of a lattice version of the 
BRST symmetry, which is currently under investigation 
ffl\ . Given such a construction and further assump- 
tions, a power-law behavior of the Green's functions is 
implied [2{|. This global BRST, of course, is not neces- 
sarily simply related to the ordinary perturbative BRST 
symmetry outside the domain of applicability of pertur- 
bation theory, and therefore statements on just the per- 
turbative BRST symmetry cannot be carried over to the 
non-perturbative domain 13 . It should be noted that in 



the Kugo-Ojima construction it is actually not relevant 
which particular non-perturbative extension of the Lan- 
dau (or linear covariant) gauge is used. The only ingre- 
dient necessary is the existence of an unbroken BRST 
charge and its algebraic properties, as can be seen di- 
rectly from the construction [7| . The existence of Gribov 
copies then implies that such a charge has necessarily 
to be constructed globally. Conversely, this implicitly 
defines a unique, non-perturbative gauge. However, an 
explicit construction would be necessary, or at least a 
proof of its existence. Until then, the existence of such 
a charge can only be assumed, as has been done in [29| . 
Whether this gauge is then identical to the absolute Lan- 
dau gauge can be checked, at least in principle, by the 
condition ([3]). An infrared singular ghost dressing func- 
tion in the absolute Landau gauge would, however, be 
already an indication of such an equivalence. 

This completes the scenario, which would provide 
a picture where both, decoupling-type solutions and 
scaling-type solutions occur, though in different gauges. 
Of course, this is only one possible scenario. 

Hence, at the current point it seems more appropriate 
to compare the lattice results in absolute Landau gauge, 
rather than in minimal Landau gauge, with those from 
functional results which exhibit a scaling behavior in the 
far infrared. This is shown in figure 1161 The agreement 
is rather good at these small volumes. 



VI. SUMMARY 

When performing gauge-fixed calculations in the non- 
perturbative regime, it is necessary to precisely specify 
how to treat the residual gauge orbit which exists due 
to Gribov-Singer copies. In particular, gauge-dependent 
correlation functions like propagators can depend on how 
the residual gauge orbit is treated. 

The results here show that this is indeed the case when 
comparing the propagators in the minimal Landau gauge 
and the absolute Landau gauge. The differences quickly 
diminish with momenta, and are quantitatively essen- 
tially irrelevant at momenta larger than 1 GeV. However, 
in the far infrared region they lead to different behav- 
iors. Nonetheless, the changes are compatible with even 
a qualitative change of the infrared properties, although 
with the currently available statistics and volumes a fi- 
nal conclusion cannot be drawn, and the effects of dis- 
cretization effects is to some extent unexplored [65J. It 
has also be demonstrated that the effects at a fixed vol- 
ume become more severe with an increasing number of 
dimensions. 

This is only the result of the comparison of two of the 
possible non-perturbative completions of the perturba- 



E. g., the symmetry may not be manifest, as Lorentz symmetry 
at finite temperature. Alternatively, the global BRST symmetry 



may be well-defined and unbroken, while the ordinary one would 
not be, similar to translation symmetry on a sphere. 
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FIG. 16: The lattice results in absolute Landau gauge from a 
64 3 lattice at (3 — 4.24 compared to results from DSE calcu- 
lations [!}. 



tive Landau gauge. As discussed, at least on the lattice 
it is possible to construct a whole one-parameter family 
of Landau gauges all corresponding to the perturbative 
Landau gauge. A possible scenario of how these may be 
related to the one-parameter family of solutions found in 
the corresponding functional calculations in the infinite- 
volume and continuum limit [29| has been presented in 
section IV Bl 

There results emphasize that in the non-perturbative 
domain it is necessary to really completely specify the 
gauge-fixing procedure, also in functional methods, as 
the different gauges can correspond to different Green's 
functions. Comparisons between different methods can 
only be made if the same gauge has been employed. 
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APPENDIX A: TECHNICAL DETAILS OF THE 
LATTICE SIMULATIONS 

The configurations have been obtained using a hybrid 
heat-bath overrelaxation procedure using a standard Wil- 
son action. Details of this procedure and the determina- 
tion of the propagators along with statistical error deter- 
mination procedures can be found in 55]. 

In total, three sets of data have been obtained. One set 
has been used in section IIIII to investigate properties of 
the adaptive stochastic overrelaxation gauge-fixing pro- 
cedure 55], which is used as an element in the method 
to fix to the absolute Landau gauge described in the 
next appendix [Bl The details of these simulations can 
be found in table Q] In this set of data, the gauge-fixing 
procedure was limited in the number of iterations, as de- 
scribed in section UTTl For each limit, an independent set 
of configurations was used, to reduce systematic errors. 

Note that for all data sets in tables [Ti l VII multiple in- 
dependent runs were always performed to reduce auto- 
correlation effects. In case of the largest volumes, usually 
less than ten configurations constitute an independent 
run. 

A second set is obtained and fixed to the minimal Lan- 
dau gauge using adaptive stochastic overrelaxation [55] ]. 
In this random gauge transformation is used to 

initialize the iterative procedure. The results from these 
configurations are used for comparison between the min- 
imal Landau gauge and the absolute Landau gauge in 
section IIVI The details of these configurations are listed 
in table [TTJ in two dimensions and in table IIIII in three 
dimensions. The number of configurations was tuned to 
obtain the exponent of the ghost propagator with a sta- 
tistical accuracy of better than 10% at the la-level. 

The third set contains the configurations which have 
been fixed to the absolute Landau gauge. Details on 
these can be found for two dimensions in table IIVI in 
table [V] for three dimensions, and in table I VII for four 
dimensions. Note that in this case it was not always pos- 
sible to achieve the desired statistical accuracy, as it took, 
for the largest volume, up to one CPU-week on a current 
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TABLE I: Number of configurations considered in the numer- 
ical simulations for the results of section Mil The value of 
the lattice spacing a has been taken from Ref. [5lj]. All cal- 
culations have been performed in three dimensions on a 40 3 
lattice at /3 = 4.2. Always 200 thermalization sweeps and 50 
sweeps between two consecutive measurements have been per- 
formed. 'Limit' is the maximum number of permitted gauge 
fixing sweeps, and 'Conf.' the number of configurations used. 
The limit oo corresponds to no restriction on the number of 
gauge-fixing sweeps, and the number in parentheses is the 
actual maximum number observed. The data in this case is 
taken from [5(| , and the first and second number of configura- 
tions refers to the number of configurations used for the ghost 
and gluon propagator, respectively. Note that for each limit 
case independent configurations have been obtained. Note 
further that in all cases multiple independent runs (about 15- 
20) have been performed to obtain the data. 
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standard processor to obtain one gauge-fixed configura- 
tion. The final statistical accuracy achieved was the best 
obtainable with the available resources. Additional infor- 
mation concerning the gauge- fixing procedure is provided 
in the tables llVllVIl The description of these can be found 
in the next appendix [B] on the gauge-fixing algorithm. 



APPENDIX B: FIXING TO THE ABSOLUTE 
LANDAU GAUGE 

Fixing to the absolute Landau gauge as specified in sec- 
tion[TT]is a hard problem. Actually it is NP-hard, similar 
to finding the ground state of a spin glass Hamiltonian, 



TABLE II: Number of configurations considered in the nu- 
merical simulations for the results in minimal Landau gauge 
of section ITVl in two dimensions. The value of the lattice spac- 
ing a has been determined analytically [79( |. To assign units 
the standard value of (440 MeV) 2 has been assigned to the 
string tension. N is the number of lattice sites in one di- 
rection of the symmetric lattice, a is the lattice spacing, L 
the physical length extent, 'Therm.' the number of thermal- 
ization sweeps, 'Sweep' the number of sweeps between two 
consecutive measurements and 'Conf.' the total number of 
configurations. Note that in all cases multiple independent 
runs have been performed to obtain the data. 
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to the best of current knowledge. Therefore, there is 
no hope to really find the absolute minimum of ((3]). In 
particular, confirming to have found the absolute min- 
imum of ([3]) is in itself a problem at least as hard as 
the original problem of finding the minimum. Therefore, 
the best attempt, which is pursued here, is to improve 
the condition (j4|). That this at least resembles the ab- 
solute Landau gauge rests on the assumption that the 
evolution of all (relevant) quantities with the value of ([3]) 
is smooth, and the limit itself is also not discontinuous. 
Hence, by getting closer and closer to the absolute mini- 
mum the propagators should become more and more like 
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TABLE III: Number of configurations considered in the nu- 
merical simulations for the results in minimal Landau gauge of 
section ITVl in three dimensions. The value of the lattice spac- 
ing a has been determined by interpolation from Ref. [5l| . 
To assign units the standard value of (440 MeV) 2 has been 
assigned to the string tension. For further details, see table 
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2.2 


5.0 


760 


76 


1078 


32 


4.24 


0.17 


1.1 


5.4 


520 


52 


899 


64 


7.09 


0.090 


2.2 


5.8 


840 


84 


1062 


40 


4.24 


0.17 


1.1 


6.8 


600 


60 


1040 


48 


4.24 


0.17 


1.1 


8.2 


680 


68 


1047 


56 


4.24 


0.17 


1.1 


9.5 


760 


76 


1053 


64 


4.24 


0.17 


1.1 


11 


840 


84 


1033 



the ones in the correct absolute Landau gauge. On the 
smallest lattices, there is a reasonable hope to have found 
the absolute minimum in most of the cases. There, no 
discontinuous behavior has been found, giving at least 
some justification for this assumption. 

It is then necessary to design an algorithm which aims 
at improving the condition (J3J). This has often been pur- 
sued using multiple restart algorithms [5^. | 5 ^[59j l, stimu- 
lated annealing [3] and other approaches [4a. l58l. I66L . 
However, the problem of gauge-fixing is, as has been dis- 
cussed in section [TTT1 highly configuration-dependent. It 
is therefore desirable to use an algorithm which is able 
to adapt to a given configuration. One possibility are 
genetic (or evolutionary) algorithms [69]. This type of 
algorithms has been used previously for gauge-fixing in 
lattice gauge theory [H, Ha,[63]- Here, a different flavor 
of such an algorithm is constructed and used throughout 
section IIVI 

The following algorithm has been performed for any 
given configuration C: 

1) Generate an initial population of Nw gauge trans- 
formations. They are obtained using an adap- 
tive stochastic overrelaxation algorithm from once 
the identity transformation and Nw — 1 random 
gauge transformations, and all fix C to the Landau 



gauge 

2) Determine the fitness of each member of the popu- 
lation by calculating ([3]) for it 

3) Remove the Nw/2 elements with the lowest fitness 
function, obtaining the parent generation V . Note 
that all digits of the fitness function are always con- 
sidered, i. e., even two members with a fitness only 
differing by machine precision fluctuations are not 
considered equal. Then construct a new generation 
by the following operations 

a) Generate 0.l5Nyy = N p new members by 
point mutation 15 . For this, select one mem- 
ber from V randomly. Then draw a random 
number R according to a Gaussian distribu- 
tion with a given width P p . Take the absolute 
value, and if the number is smaller than 1 set 
it to one. Replace at R randomly chosen lat- 
tice points x the local gauge transformation 
matrix g{x) by a randomly chosen SU(2) ma- 
trix. Permit multiple replacements. However, 
P p should be chosen sufficiently small. Here, 
it was selected to be 1% of the number of lat- 
tice sites. Use the so generated initial gauge 
transformation as seeds to obtain once again 
a gauge transformation of C to Landau gauge 

b) Generate 0.2iVyi/ = N c new members by cross- 
ing. For this purpose select two different mem- 
bers from V randomly. Then construct a new 
member by mixing these two parent members. 
For this, select at each lattice site randomly 
from which parent to take the matrix g(x). 
Use this new member as a seed to obtain a 
gauge transformation of C to Landau gauge 

c) Add an additional number N m of new mem- 
bers such that the population again consists of 
Nw members. These are again randomly ini- 
tialized gauge transformations, which are then 
used to construct gauge transformations of C 
to Landau gauge 

4) Determine the fitness of all members according to 
([3]). If none of the new members are better than 
those of the set V, terminate the algorithm and use 
the fittest member for the gauge transformation to 
absolute Landau gauge. If a better new member ex- 
ists, go to point 3 and repeat until the termination 
criterion is satisfied. A new member is also consid- 
ered better if it only differs from the best so far in 
fitness by an amount of the order of the machine 
precision 



Landau gauge refers here to the (local) minimum of (|3j which is 
reached by using adaptive stochastic overrelaxation [551 ■ 
Non-integer numbers are always rounded mathematically, but set 
to one if zero. 
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TABLE IV: Number of configurations considered in the numerical simulations for the results in absolute Landau gauge of 
section HVl in two dimensions. For details see table HT1 The additional parameters N m , N p , P p , N c , N m are explained in 
appendix [B] N g is the average number of generations, N x the maximum number of generations, and iV m i n , Nf, and JV max are 
the minimum number, the average number per configuration, and the maximally needed number of gauge- fixings, respectively. 



N 


a 


a [fm] 


a- 1 [GeV] 


L [fm] 


Therm. 


Sweep 


Conf. 


N w 


N p 


Pv 




N m 


N g 


N x 


iVmin 


N f 


N max 


4 


38.7 


0.089 


2.2 


0.36 


140 


14 


1027 


8 


1 


0.16 


1 


2 


2.36(2) 


6 


12 


13.44(8) 


28 


4 


10 


0.18 


1.1 


0.72 


140 


14 


1027 


8 


1 


0.16 


1 


2 


2.76(3) 


6 


12 


15.0(1) 


28 


10 


38.7 


0.089 


2.2 


0.89 


200 


20 


1208 


10 


1 


1.0 


2 


2 


2.55(2) 


9 


15 


17.8(1) 


50 


16 


38.7 


0.089 


2.2 


1.4 


260 


26 


1049 


12 


1 


2.6 


2 


3 


2.58(2) 


7 


18 


21.5(1) 


48 


10 


10 


0.18 


1.1 


1.8 


200 


20 


1027 


10 


1 


1.0 


2 


2 


2.65(3) 


8 


15 


18.3(2) 


45 


20 


38.7 


0.089 


2.2 


1.8 


300 


30 


1041 


14 


2 


4.0 


2 


3 


2.75(3) 


7 


21 


26.3(2) 


56 


26 


38.7 


0.089 


2.2 


2.3 


360 


36 


1130 


16 


2 


6.8 


3 


3 


2.53(2) 


8 


24 


28.2(2) 


72 


16 


10 


0.18 


1.1 


2.9 


260 


26 


1037 


12 


1 


2.6 


2 


3 


2.53(2) 


6 


18 


21.2(2) 


42 


34 


38.7 


0.089 


2.2 


3.0 


440 


44 


1032 


18 


2 


12 


3 


4 


2.71(3) 


7 


27 


33.4(8) 


70 


20 


10 


0.18 


1.1 


3.6 


300 


30 


1119 


14 


2 


4.0 


2 


3 


2.63(2) 


7 


21 


25.4(2) 


56 


44 


38.7 


0.089 


2.2 


3.9 


540 


54 


1061 


20 


3 


19 


4 


3 


2.59(3) 


7 


30 


35.9(3) 


80 


26 


10 


0.18 


1.1 


4.7 


360 


36 


1065 


16 


2 


6.8 


3 


3 


2.59(3) 


6 


24 


28.7(2) 


56 


56 


38.7 


0.089 


2.2 


5.0 


660 


66 


1045 


24 


3 


31 


4 


5 


2.63(3) 


7 


36 


43.6(4) 


96 


34 


10 


0.18 


1.1 


6.1 


440 


44 


1032 


18 


2 


12 


3 


4 


2.55(3) 


7 


27 


32.0(3) 


70 


68 


38.7 


0.089 


2.2 


6.1 


680 


68 


1130 


28 


4 


46 


5 


5 


2.47(2) 


6 


42 


48.6(3) 


98 


80 


38.7 


0.089 


2.2 


7.1 


900 


90 


1084 


32 


4 


64 


6 


6 


2.52(2) 


6 


48 


56.3(3) 


112 


44 


10 


0.18 


1.1 


7.9 


540 


54 


1006 


20 


3 


19 


4 


3 


2.46(2) 


6 


30 


34.6(2) 


70 


104 


38.7 


0.089 


2.2 


9.3 


1140 


114 


2078 


36 


5 


108 


7 


6 


2.50(2) 


7 


54 


63.0(4) 


144 


56 


10 


0.18 


1.1 


10 


660 


66 


1083 


24 


3 


31 


4 


5 


2.56(2) 


7 


36 


42.7(2) 


96 


112 


38.7 


0.089 


2.2 


10 


1220 


122 


1068 


40 


6 


125 


8 


6 


2.50(3) 


7 


60 


70.0(6) 


160 


68 


10 


0.18 


1.1 


12 


780 


78 


1014 


28 


4 


46 


5 


5 


2.47(2) 


6 


42 


48.6(3) 


98 


136 


38.7 


0.089 


2.2 


12 


1460 


146 


1028 


46 


6 


185 


9 


8 


2.49(2) 


7 


69 


80.3(5) 


184 


150 


38.7 


0.089 


2.2 


13 


1600 


160 


459 


52 


7 


225 


10 


9 


2.49(3) 


6 


78 


90.7(8) 


182 


80 


10 


0.18 


1.1 


14 


900 


90 


2136 


32 


4 


64 


6 


6 


2.45(2) 


8 


48 


55.2(3) 


144 


166 


38.7 


0.089 


2.2 


15 


1760 


176 


326 


58 


8 


276 


11 


10 


2.52(5) 


6 


87 


102(1) 


203 


184 


38.7 


0.089 


2.2 


16 


1940 


194 


66 


64 


9 


339 


12 


11 


2.45(9) 


5 


96 


110(3) 


192 


104 


10 


0.18 


1.1 


19 


1140 


114 


1096 


36 


5 


108 


7 


6 


2.44(2) 


6 


54 


61.9(4) 


126 


112 


10 


0.18 


1.1 


20 


1220 


122 


2105 


40 


6 


125 


8 


6 


2.45(2) 


10 


60 


69.0(4) 


220 


136 


10 


0.18 


1.1 


24 


1460 


146 


1414 


46 


6 


185 


9 


8 


2.46(2) 


7 


69 


79.6(5) 


184 


150 


10 


0.18 


1.1 


27 


1600 


160 


748 


52 


7 


225 


10 


9 


2.51(3) 


7 


78 


91.3(8) 


208 


166 


10 


0.18 


1.1 


30 


1760 


176 


267 


58 


8 


276 


11 


10 


2.52(6) 


6 


87 


102(2) 


203 


184 


10 


0.18 


1.1 


33 


1940 


194 


53 


64 


9 


339 


12 


11 


2.36(8) 


4 


96 


108(3) 


160 



Note that this algorithm does not guarantee to termi- 
nate after a number of generations. However, in practice 
it always terminates. In most cases, already the second 
generation did not yield any improvement, but the max- 
imum number of generations N x has been as large as 
10 for the systems investigated here. The average num- 
ber of generations N g is always between two and a bit 
more than three, yielding therefore a total number of per- 
formed gauge- fixings between 3/2Nw and roughly 2Nw 
Given that the number of Gribov copies grows exponen- 
tially with the number of lattice sites, N\y also has to 
grow with the number of lattice sites. Strict exponen- 
tial growth is at this time prohibited by computational 
resources. The values chosen in practice are listed in ta- 



bles HvED 

Note that the numbers N\y, N p , N c , N m , and P p are 
far from optimized, and result only from experiments on 
small lattices. This algorithm can and should be im- 
proved, if fixing to the absolute Landau gauge turns out 
to be important for applications. Furthermore, the struc- 
ture of the algorithm is itself subject to optimization, and 
may possibly be more efficient if other structures were ad- 
mitted as well [69[ . Finally, there exists some knowledge 
on the structure of Gribov copies (8l| . In particular, they 
exhibit a domain-like structure. Incorporating this prior 
knowledge in, e. g., the crossing of two gauge transfor- 
mations may furthermore improve the efficiency of the 
algorithm. 



[1] R. J. Rivers, "Path integral methods in quantum field bridge monographs on mathematical physics). 

theory," Cambridge, UK: Univ. Pr. (1987) 339 p. (Cam- [2] M. Bohm, A. Denner and H. Joos, "Gauge Theories 
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TABLE V: Number of configurations considered in the numerical simulations for the results in absolute Landau gauge of section 
IIVI in three dimensions. See table [TTT1 and HVl as well as appendix [Bl for details. 



N 


13 


a [fm] 


a' 1 [GcV] 


L [fm] 


Therm. 


Sweep 


Conf. 


N w 


N p 


P P 


N c 


N m 


N g 


N x 


iVmin 


N f 


N: max 


4 


7.09 


0.090 


2.2 


0.36 


240 


24 


1092 


8 


1 


0.64 


1 


2 


2.70(2) 


6 


12 


14.8(2) 


28 


4 


4.24 


0.17 


1.1 


0.68 


240 


24 


1270 


8 


1 


0.64 


1 


2 


2.67(2) 


7 


12 


14.7(2) 


32 


8 


7.09 


0.090 


2.2 


0.72 


280 


28 


1115 


10 


1 


5.1 


2 


2 


2.69(3) 


9 


15 


18.5(2) 


50 


12 


7.09 


0.090 


2.2 


1.1 


320 


32 


1115 


12 


1 


17 


2 


3 


2.68(3) 


7 


18 


22.1(2) 


48 


8 


4.24 


0.17 


1.1 


1.4 


280 


28 


1196 


10 


1 


5.1 


2 


2 


2.73(3) 


7 


15 


18.7(2) 


40 


16 


7.09 


0.090 


2.2 


1.4 


360 


36 


1106 


14 


2 


41 


2 


3 


2.62(3) 


8 


21 


25.3(2) 


63 


20 


7.09 


0.090 


2.2 


1.8 


400 


40 


1144 


16 


2 


80 


3 


3 


2.64(2) 


8 


24 


29.1(2) 


72 


12 


4.24 


0.17 


1.1 


2.0 


320 


32 


1109 


12 


1 


17 


2 


3 


2.54(2) 


6 


18 


21.2(1) 


42 


24 


7.09 


0.090 


2.2 


2.2 


440 


44 


1093 


18 


2 


138 


3 


4 


2.56(3) 


7 


27 


32.0(3) 


72 


16 


4.24 


0.17 


1.1 


2.7 


360 


36 


1041 


14 


2 


41 


2 


3 


2.54(3) 


8 


21 


24.8(2) 


63 


32 


7.09 


0.090 


2.2 


2.9 


520 


52 


1025 


20 


3 


328 


4 


3 


2.68(3) 


8 


30 


36.8(3) 


90 


20 


4.24 


0.17 


1.1 


3.4 


400 


40 


1302 


16 


2 


80 


3 


3 


2.60(3) 


8 


24 


28.8(2) 


72 


40 


7.09 


0.090 


2.2 


3.6 


600 


60 


1005 


24 


3 


640 


4 


5 


2.73(3) 


8 


36 


44.8(4) 


108 


24 


4.24 


0.17 


1.1 


4.1 


440 


44 


1095 


18 


2 


138 


3 


4 


2.57(3) 


7 


27 


32.1(3) 


72 


48 


7.09 


0.090 


2.2 


4.3 


680 


68 


906 


28 


4 


1106 


5 


5 


2.84(3) 


8 


42 


53.8(4) 


126 


56 


7.09 


0.090 


2.2 


5.0 


760 


76 


657 


32 


4 


1756 


6 


6 


2.83(4) 


8 


48 


61.3(6) 


144 


32 


4.24 


0.17 


1.1 


5.4 


520 


52 


1072 


20 


3 


328 


4 


3 


2.71(3) 


9 


30 


37.1(3) 


100 


64 


7.09 


0.090 


2.2 


5.8 


840 


84 


167 


36 


5 


2621 


7 


6 


3.02(7) 


6 


54 


72(1) 


126 


40 


4.24 


0.17 


1.1 


6.8 


600 


60 


1026 


24 


3 


640 


4 


5 


2.72(3) 


7 


36 


44.6(4) 


96 


48 


4.24 


0.17 


1.1 


8.2 


680 


68 


929 


28 


4 


1106 


5 


5 


2.86(3) 


8 


42 


54.0(4) 


126 


56 


4.24 


0.17 


1.1 


9.5 


760 


76 


522 


32 


4 


1756 


6 


6 


2.87(4) 


7 


48 


61.9(6) 


128 


64 


4.24 


0.17 


1.1 


11 


840 


84 


291 


36 


5 


2621 


7 


6 


3.07(6) 


7 


54 


73(1) 


144 



TABLE VI: Number of configurations considered in the numerical simulations for the results in absolute Landau gauge of 
section HVl in four dimensions. See table 11111 and ITVl as well as appendix [Bl for details. The scale has been fixed according to 
0. 



N 




a [fm] 


a" 1 [GcV] 


L [fm] 


Therm. 


Sweep 


Conf. 


N w 


N p 


Pv 


N c 


N m 


N g 


N x 


iVmin 


Nf 


N: max 


4 


2.45 


0.099 


2.0 


0.40 


340 


34 


1020 


8 


1 


2.6 


1 


2 


2.63(3) 


7 


12 


14.5(1) 


32 


6 


2.45 


0.099 


2.0 


0.59 


360 


36 


1191 


10 


1 


13 


2 


2 


2.60(3) 


7 


15 


18.0(2) 


40 


8 


2.45 


0.099 


2.0 


0.79 


380 


38 


1013 


12 


1 


41 


2 


3 


2.59(2) 


8 


18 


21.5(1) 


56 


4 


2.2 


0.21 


0.94 


0.84 


340 


34 


1132 


8 


1 


2.6 


1 


2 


2.63(2) 


7 


12 


14.5(1) 


32 


12 


2.45 


0.099 


2.0 


1.2 


420 


42 


1091 


14 


2 


207 


2 


3 


2.64(3) 


9 


21 


25.5(2) 


70 


6 


2.2 


0.21 


0.94 


1.3 


360 


36 


1191 


10 


1 


13 


2 


2 


2.57(3) 


8 


15 


17.9(2) 


45 


16 


2.45 


0.099 


2.0 


1.6 


460 


46 


652 


16 


2 


655 


3 


3 


2.65(4) 


8 


24 


29.2(3) 


72 


8 


2.2 


0.21 


0.94 


1.7 


380 


38 


1013 


12 


1 


41 


2 


3 


2.56(3) 


6 


18 


21.4(2) 


56 


20 


2.45 


0.099 


2.0 


2.0 


500 


50 


148 


18 


2 


1600 


3 


4 


2.72(9) 


10 


27 


33.5(7) 


100 


12 


2.2 


0.099 


0.94 


2.5 


420 


42 


1081 


14 


2 


207 


2 


3 


2.65(3) 


7 


21 


25.6(2) 


56 


16 


2.2 


0.21 


0.94 


3.4 


460 


46 


277 


16 


2 


655 


3 


3 


2.78(6) 


8 


24 


30.2(5) 


72 


20 


2.2 


0.21 


0.94 


4.2 


500 


50 


59 


18 


2 


1600 


3 


4 


2.7(1) 


5 


27 


33.3(8) 


54 
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